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Abstract 


The equations of motion of a geodetic satellite in the earth's gravi- 
tational field expressed by gravity anomalies require the evaluation, amongst 
others, of the partial derivatives of the disturbing force with respect to 
individual gravity anomalies. These derivatives would be in error if evaluated 
using coordinates at the center points of the mean gravity anomaly blocks. 

This report discusses how these blocks should be subdivided so that the 
partial derivatives could be numerically evaluated for each subdivision, and 
then finally meaned to give the value representative of the whole blocks, with 
accuracies better than 2-3% for all blocks. The number of subdivisions is 
large for the blocks nearest to the satellite subpoint and decreases away from 
it. The actual values of this spherical distance and the actual subdivision of 
the mean gravity anomaly blocks has been determined numerically for 184 
15° x 15° equal area blocks. Satellite heights above the earth of 400 km, 800 km 
and 1600 km have been considered. The computer times for the suggested 
scheme have been compared with alternative solutions. 
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1. Introduction 


The predominant force acting on artificial earth satellites used 
for geodetic purposes is that of the earth's gravity field. If we choose to 
represent it by a set of mean gravity anomalies over specified blocks and 
referred to a defined reference surface, we need to evaluate the effect of 
each gravity anomaly at the given satellite position. As this has to be 
repeated for each satellite position being considered, for example in the 
numerical integration approach of orbital and trajectory analysis, any 
simplification in the practical evaluation consistent with the required 
accuracy would be significant. 

The following specific problems have been examined in this paper: 

(a) Is it adequate to consider the effect from the center of the block 
over which the mean gravity anomaly is given, or should the 
effect be meaned over several points in the block? If so, how 
should these points be chosen? 

(b) Is it adequate to have the computations being meaned over points 
in the block for only a few blocks near the satellite subpoint; and 
if so, up to what distance from the satellite subpoint? 

(c) Will it be possible to ignore the effect of some blocks altogether, 
which are far removed from the satellite subpoint? 

(d) How will the above conclusions vary with change in the height of 
the satellite? 

2. The Basic Equations 


Following Rapp (1971), the equations of motion of the satellite in an 
inertial coordinate system (x, y,z) at time t, measured from an initial epoch 
to, in terms of acceleration components x,y,z, may be expressed in a general 
form as: 



• • • 

where Xo ,y 0 , z 0 and Xo,y 0 »z 0 are the initial position and velocity components 
of the satellite at epoch t 0 , and Agi, Ag s , . . . Ag n are the mean gravity 
anomalies. 

Using {3 k for any one of the individual gravity anomalies or the initial 
position and velocity components (x 0 ,y 0 ,z 0 , Xo,y 0 ,z 0 ), the variational equations 
with respect to fa may be expressed as: 
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We will confine our discussion to the evaluation of partial derivatives 
with respect to the gravity anomalies and use the notation:’ 




3h 


(3) 


where the evaluation of C x , C y ,C Z will be done for each individual gravity 
anomaly Ag x , A g 3 , . . . Ag n . 

3. Equations for C T , C v , C r 


Considering the equations of motion (1) of the satellite only due to the 
earth’s gravity potential composed of the normal part U and the disturbing part 
T, we may write: , ... 


x = 

ax 


au 

ax 


a t 
ax 


(4) 


and similar expressions for y , z , which when combined with equations (3) and 


2 



(2) give: 



and similar expressions for C y ^ ,, C z for an individual gravity anomaly Ag t . 

As the partial derivatives of the disturbing potential with respect to 
the satellite position in the x, y, z coordinate system may be expressed with 
respect tor , i|> , and X.,' being respectively the geocentric radius vector, 
latitude and longitude of the satellite; e.g. , 
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and as the value of 

a r 

, a* , 

ax 

may be taken from literature (eg. 


Heiskanen and Moritz, 1967) in terms of the generalized Stokes' function, 
we may express the value of C x , C y , C z following Rapp (1971) as: 



S2L dX_ 3X_ ax 

where $ x , Sy , 5z have been used in place of $x , 3y , Sz these 
being equal at any given epoch as a, the Right Ascension of the satellite is 
given in terms of X and 0 , the Greenwich Side rial time by 


O' = X + 9 


(8) 
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The values of the coefficients A 4 , B t , Ci corresponding to a given 
gravity anomaly Ag t may be expressed as: 


A t 


LtJlL da 
4n cir 
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tR cos 
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dS(r-, »*) 
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(9) 


where R is the radius of a sphere approximating the earth, da is the. area of 
the gravity anomaly block/sub block \|f * is the spherical distance and a* is 
the azimuth from the satellite subpoint to the gravity anomaly block/subblock 
and S( r ,i|r*) is the generalized Stokes’ function. Its value and that of its 
partial derivatives with respect to rand ^ * is given as in Heiskanen and 
Moritz (1967) by: 


S( r *) 


as 

9 r 
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as 
a* * 


= |-3 + |- -32m L-tcosf!+D . 

L D D ° Dsin a i|r* 2 


R 


t 2 )* 


where t = — and D = (l-2t cos i|t* + t 2 ) 
r 

The partial derivatives of r, a, f with respect to the satellite position 
are obtained from: 


x = rcos i|rcos oe 
y = rcos \(f sin a 
z = r sin i|r 
r 3 = x 3 + y 3 + z a 


( 11 ) 
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Finally, the values of i|r * and a* may be obtained in terms of the, geocentric 
latitude i|r and longitude X of the satellite subpoint and the corresponding values 
i|r g and X g of the gravity anomaly block/sub block by: 


cos i|r * = sin i|r sin i|r g + cos i|r cos t|r g cos ( \ - X) 

sin a* = -£2g jL sin (X K -X) 

sin i|r * , (12) 

cos a* = cos sin i|t g - sin f cos cos(X e - X) 

sin f * 


The value of C Xj , C y , C z for each of the gravity anomalies Ag t 
(Agi, Ag s , . . . Ag n ) for a given satellite position (x,y,z) at epoch t is thus 
given by substituting equations (9) to (12) in equation (7). 

We have confined our discussion in this paper to the evaluation of 
C x , Cy , C z numerically. The purpose of the investigation has been to compute 
them accurately enough in an optimum manner, thereby saving computer time, 
by omitting excessive computations, which do not lead to appreciable increase 
in accuracy. The computations were done on an IBM System/370 Model 165 
computer. 

4. Test Data 


The gravity anomalies used to define the earth's gravity field were 
184 15° x IE? equal area mean anomalies, as listed by Obenson [1970, pp. 

127 - 128] and used by Haverland [1971], 

The satellite orbit was generated using a slightly modified version of 
a Cowell orbit generation program using an eleventh order Cowell integration 
witha fixed step size of 60 seconds. The earth's gravitational field alone was 
considered using potential coefficients C 3 0 and C 4 0 , as in Rapp (1971). 

The computer program as developed by Haverland (1971) and subsequently 
revised by Rapp (unpublished) was used. A few modifications were made but the 
orbital parameters were retained, except for different values for the semi-major 
axis to give approximate satellite heights above the earth of 400 km, 800 km, 
and 1600 km to cover the usual range of geodetic satellites. For still higher 
satellites, the conclusions drawn in this paper will be on the safe side. 
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The other orbital parameters were: 


eccentricity 
inclination 
argument of perigee 
longitude of ascending node 
mean anomaly at epoch 


e = 0.07061... 
i = 59°. 386... 

Co = 312° . 74. . . 

Q = 263°. 74. . . 

M c = 0.23176... revolutions 


5. Block Subdivisions 


The most direct evaluation of C* , C y ,C t would take place by evaluating 
equation (9) using coordinates ofthe center of the 15° block only. Such an 
evaluation may suffer from a numerical integration error. Consequently, we 
first examine a subdivision of the block considering 2 sub blocks and 4 sub 
blocks respectively. The value of C x ,C y ,C Z representative of the whole block 
would be obtained as a mean of the sub block values. These means could then 
be termed as 2 point mean and 4 point mean respectively, and could be used 
in place of one point (center point) value of the whole IS 3 x Iff block. As the 
latitudinal extent, Acp of the blocks was uniformly Iff, but as the longitudinal 
extent AX , increased towards the poles up to 120° from the 15° value at the 
equator to keep the blocks nominally as equal area, the scheme of subdivision 
of the blocks into sub blocks was as below: 



The value of C* , C y , C 2 was computed based oh a 1 point (center ) 
evaluation and as 2 point mean and 4 point mean values for all the 184 
15° x 15° blocks. The differences (2 point mean - 1 point value) and ( 4 point 
mean - 1 point value) were computed individually for each mean gravity anomaly 
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block. A root mean square value was then obtained for these differences over 
all the 184 blocks and compared against the root mean square value over the 
184 blocks of the 2 point mean and 4 point mean values of C„ ,C y , C z , and 
expressed as a percentage. 

The results are given in Table 1 for four different satellite heights. 

It was found that the 2 point mean value, when compared with the 4 point mean 
value, was no better than the one point value, but worse in most cases. This 
is perhaps due to the assymetric position of the centers of the sub blocks 
in the case of 2 point means, as compared to the center of the whole block. 
Accordingly, 2 point means have not been considered further. For the same 
reasons, no other assymetric subdivisions of the block was considered. 


Following these first results, further symmetric subdivisions were 
chosen to obtain more refined values. Specifically a 9 point and a 16 point 
mean was evaluated. This subdivision is shown below: 
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The results in Table 1 show the effect of increasing the number of 
sub blocks on the values of C x , C y , C * . It is reasonable to expect that the 
16 point means are more accurate that 9 point means, which are more accurate 
than 4 point means* etc. However, as root mean square values over all the 
184 blocks have been considered in Table 1, the root mean square value of 
differences (16 point mean - 1 point value), (16 point mean - 4 point value), etc. 
come out to be rather large. As we see later, these differences are predominantly 
large for individual gravity anomaly blocks nearest to the satellite subpoint, 
and are much smaller away from it. 


The results in Table 1 may be used as a guide in estimating the error 
in computing C xf C y ,C z when using a 1,4, or 9 point mean if we accept the 16 
point mean as being correct. 
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Effect of Number of Block Subdivisions on Accuracy of Computation 
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6. Number of Blocks for Computation of 16 Point/ 9 Point Means 


The predominant difference between the 16, 9 and 4 point means occur 
for those mean gravity anomaly blocks, which are nearest to the satellite sub- 
point, within a spherical distance ty* (equation (12)) of 10° - 20° . 

Tables 2, 3 and 4 show respectively for satellite heights above earth of 
about 400, 800 and 1600 km, the values of C* , C y , C z for 16,9,4 point means and 
center (1 point) value of all mean gravity anomaly blocks within a spherical 
distance i|r* of 20° from the satellite subpoint. 

It is clear from these results that center point value cannot be used 
for i|f*< 20, and even 4 point mean is in considerable error for ijKlS 0 or so. 

For the first nearest block with i|r* <5°, even the 9 point mean appears to be 
in error by 1% to 4%. 
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Values for Nearest Blocks ^*<20°Satellite Ht. e®400 km 
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Values for Nearest Blocks < 20° Satellite Ht. 800 km, 
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We may therefore conclude that for errors due to subdivision of 
blocks not to exceed about 2%, a 16 point mean should be taken when the 
center (cp g , X g ) of mean gravity anomaly block is less than 10° from the 
satellite subpoint, and' a 9 point mean when 10° < ♦*< 15° . For determining 
the value of up to which 4 point mean should be taken, we should examine 
further results and this has been done in Section 7. 

We also find from Tables 2, 3, 4 that values of , Cy , C g decrease 
with height of the satellite and the percentage value . of difference of means 
(qg. 16 pt. mean - 4 pt.mean) also becomes less. The limits of <10° and 
10° - Iff for 16/9 point means however appear to hold for satellite heights 
from 400 km to 1600 km. 

Apart from the above strictly numerical point of view, we may also 
consider the optimum choice of spherical radius i|r * for 16pt/9pt mean 
computation from the symmetry and the number of mean gravity anomaly 
blocks falling within the chosen value of i|i*. In Figure 1, we have considered 
three cases for the location of the satellite subpoint in relation to the mean 
gravity anomaly blocks; 

(1) At the common corner of 4 blocks, 

(2) In the middle of the common edge of 2 blocks, 

(3) hi the center of a block. 

Further, as the longitudinal extent of the blocks increase away from equator, 
the equatorial case above has been shown. At other locations of the satellite 
subpoint away from the equator, the number of mean gravity anomaly blocks 
with their centers (\|f g , X g ) within the specified i|f* will only be less. 

For each of the three cases, circles of radii i|f* = 5° to 20° after 
every 2.°5, have been drawn, and the number of mean gravity anomaly blocks, 
whose centers fall within these are shown in Table 5. 

TABLE 5 

No. of Mean Gravity Anomaly Blocks for given iff*(5°ojf*£20 0 ) 


ijr * 

Case 1 • 

Case 2 

Case 3 

7.5° - 10° 

0 

2 

1 

12. 5° 

4 _ ? .. 

2 

1 

17.5° 

4 

6 

5 

to 

o 

o 

4 

6 

5 
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Figure 1 From 5° To 20° From Satellite Subpoi 


4 


CASE: I. 

SATELLITE SUBPOINT 
AT COMMON CORNER 
OF 4 BLOCKS 


CASE 2. 

SATELLITE SUBPOINT 
AT MIDPOINT COMMON 
EDGE OF 2 BLOCKS 


CASE 3. 

SATELLITE SUBPOINT 
AT CENTER OF A BLOCK 




From this, and in view of the numerical. results in Table 1, we 
may choose to have 16 point mean for 0°<i|f *<10° and 9 point mean 10°^ 
t ^ 20° , the slight increase in ^ * not appreciably increasing the number) 
of blocks. .... • . !-• . 

The maximum number of mean gravity anomaly blocks requiring 
computation of 16 pt/9 pt mean would then be 1 or 2 . and an additional 2 to 
4 respectively, depending on the location of the satellite subpoint. When 
this happens to be within 5° of the pole, there will be only 3 blocks for 
16/9 pt mean computation. 

7. Number of Blocks for Computation of 4 Point Means 


We now examine in Tables 6, 7 and 8 for satellite heights of approx- 
imately 400, 800 and 1600 km, the values of C* , C y , C t for those mean 
gravity anomaly blocks, whose centers (c%,X g ) are at spherical distance 
of 20° - 45° from the satellite subpoint. For satellite height of about 400 km, 
1 point value agrees with 16 point mean within about 2. 5%, only after <|f*>40°. 
For the same tolerance, f* >35° for satellite height of about 800 km, and 
♦ * >30° for satellite height of about 1600 km. We may therefore compute 
4 point means for the mean gravity anomaly blocks whose centers are within 
the above spherical radius from the satellite subpoint and compute only 
the center (1 point) value after that. : 
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Values for 20 < t* < 45° Satellite Ht. « 400 km 
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Values for 20° < ♦* < 45° Satellite Ht. « 1600 km 
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We may now again consider the optimum choice of ijr * from the 
symmetry and number of blocks falling within a specified value of i|f *. 
Figures 2,3,4 show circles with radii i|f * = 30°, 35°, 40°, 45° for the 
cases 1,2,3 of the satellite subpoint location, as discussed in Section 
6. These figures are for the satellite subpoint being on the equator. 

Figure 5 shows the location of mean anomaly gravity blocks around either 
pole. The number of blocks when the satellite subpoint is in mid-latitudes 
will be somewherebetween the equatorial and polar cases. The results for 
the latter two cases have been summarized in Table 9 below. 


TABLE 9 

Number of Mean Gravity Anomaly Blocks forgiven \|r* 1 30°^i|f *^45°) 


jjf * 

Fig. 2 

Fig, 3 

Fig. 4 

Fig. 5 

°o 

CO 

12 

12 

13 

12 

35° 

16 

16 

18 

12 

40° 

24 

22 

21 

27 

4^ 

24 

26 

24 

27 
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Figure 2 


30° ^ i|r * s 45° Equatorial Location Case 1 
Satellite Subpoint at Common Corner of 4 Blocks 
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Figure 3 30° * f * s 45 ° Equatorial Location Case 2 

Satellite Subpoint at Midpoint of Common Edge of 2 Blocks 

21 



Figure 4 


30° ^ ♦ * s 45 0 


Equatorial Location 


Case 3 


Satellite Subpoint at the Center of a Block 
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Clearly, for satellites, with heights above 1600 km or so, it will 
be adequate to take 4 point means for iJj* £ 30°, the number of blocks being 
12 to 13. This number includes the 4 to 6 blocks for which 9/16 point means 
will be taken as discussed in Section 6. For satellites of height of about 800 
km, we may take 4 point means for p £ 35°, the number of blocks being 12 
to 18, usually 16. However, for satellites of still lower heights, the number 
of blocks for p ^ 40° is 21 to 27, which is not significantly different from the 
number of blocks for p ^45°, being 24 to 27. 

It is therefore worthwhile to compute 4 point means for all blocks, 
whose centers are within a spherical radius 20° < P^ 45° from the satellite 
subpoint, if the height of the satellite is lower than 800 km. For satellite 
heights 800 km to 1600 km, 20°<P ^35°. For higher satellites, 4 point 
means may be taken for 20° < t|t * ^ 30°. 

8. Computation of Blocks at >135° 


We now examine if the values of C x , C y , C z for the mean gravity anomaly 
blocks, which are far away from the satellite subpoint, say at spherical distance 
P >13^, are so small as compared to the average values, that they may be 
ignored, i.e. not computed, and assumed to be zero. 

The root mean square value of C x , C y ,C Z was accordingly first examined 
for all the 184 blocks, and then for the remaining blocks after respectively 
neglecting the blocks with p >16Jf*, >150°, > 135°. The results for satellite 
height of about 800 km are given in Table 10 below. 


TABLE 10 


Variation of RMS Value of C„ , C y , C z as Blocks with i|r>135° are Neglected 

Satellite Height « 800 km 


Root M 

ean Square Values 

Percentage Variation in RMS Values 


Pup to 
180° 

pup to 
165° 

l|l *up to 
150° 

'Pup to 
13 E? 

Pup to 165° 



Remarks 


gag 

.0259 


0.7 

3.2 

7.6 

rc±7187km 
h~ 816 km 

No. of Blocks ] 

Neglecte 

d 

3 

13 

27 

t = 1.0 min 




0.4 

2.4 

5.9 


No. of 

Blocks] 

Neglecte 

d 

3 

13 

27 





0.6 

2.7 

6. 5 


No. of 

Blocks 

Neglecte 

id 

3 

* 13 

27 
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We therefore find that the values of C„, 0 y , are not negligibly small 
for large values of i|f*» when compared with their average values. This is 
seen more explicitly by listing the values of C x , C y ,C Z in Table 11 for the 
blocks at spherical distance i(i* exceeding 135°. We find that the values do 
not decrease uniformly with increase in $ *, and there are several values, which 
are more than half of the root mean square . value. 

Figure 6 displays these results graphically, and in particular, shows 
that the values are in fact larger than for other blocks, for example, for blocks 
with t * from 75° to 90° , whose values have also been shown in the same graph. 
The results for other satellite heights are similar and have thus not been shown. 

We thus cannot dispense with the computation of C*, C y , C t values for 
large values of i|< *, though it is adequate to compute them only for the center 
(1 point value) of the block, as discussed in Section 7. 
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TABLE 11 

v 

c x , Cy, C z Values for t* >135° 
Satellite Ht. «, 800 km 


Block 

No. 

Sph. 

Dist. 

i|i*Deg. 

B 

■ 

C z 

Block 

No. 

Sph. 

Dist. 

'I'* Deg 

/ C x 

• 

C y 

C z 

Remarks 

10 

135. 2 

-.0004 

.0195 

■■ 

19 

150.9 

-.0154 

RS9 

-.0085 

r =-7187 km 

18 

136.4 

-.0160 

.0101 


8 

151.0 

-.0113 

mm 

-.0061 

h = 816 km 

2 

136.9 

-.0088 

.0187 


60 

151. 6 

.0031 

IBBI 

-. 0180 

t = 1.0 min 

23 

139. 1 

.0031 

.0193 


22 

153.6 

-.0015 

.0222 

-. 0094 


7 

139. 6 

-.0139 

.0157 


36 

154.1 

-.0161 

.0131 

-. 0140 


35 

139.8 

-.0175 

.0079 


57 

154.4 

-. 0128 

.0101 

-. 0192 


61 

139.8 

.0081 

.0146 


59 

160.7 

-.0026 

.0i58 

Bml 


. 83 

141. 7 

.0050 

.0091 

-.0183 

39 

161.4 

-.0011 

.0211 

MBM 


80 

142.4 

-.0096 

.0031 

-.0187 

58 

162.1 

-.0082 

.0137 

-.0218 


56 

143.0 

-.0158 

.0056 


20 

163.9 

-.0121 

.0120 

-.0125 


40 

146.9 

.0042 

.0196 


21 

165. 5 

-.0070 

.0224 

-.0129 


82 

147.0 

-.0003 

. 0084 


37 

168.7 

-.0123 

.0175 

-.0180 


81 

9 

147.3 

148.7 

-.0054 

-.0056 

.0063 

.0227 


38 

175. 7 
RMS 
Value 

-.0070 

.0251 



.0203 

.0293 

-.0189 
' .0302 
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Satellite Ht. 800 km 


9. Timing of Computer Runs 


Computer runs were made on IBM System/370 Model 165 computer 
for the computation of C x ,C y ,C z for a specified position of the satellite for 
all the 184 15° x 15° mean gravity anomaly blocks as per equations (9) to 
(12) substituted in equation (7). The actual CPU time was noted for the 
following cases* 

(a) Center point value of C x ,C y ,C z for all 184 blocks 

(b) Four point mean value of C x ,C y , C z for all 184 blocks 

'c) Nine point mean value of C x , C y ,C z for all 184 blocks 

(d) Sixteen point mean value of C x , C y , C z for all 184 blocks. 

The CPU execution time was actually obtained for 10 loops each of 
(a) to (d) above, and after dividing by 10, the time was respectively 0.034, 
0.155, 0. 334,and 0.597 seconds. These times do not include the time for the 
preliminary/subsequent computation, or for any write statements etc. , which 
are the same for any of the cases (a) to (d). The difference in time is solely 

due to the computation being done for 1,4,9 or 16 points in the blocks, and 

obtaining a mean value. 

The CPU execution time was then similarly obtained for computing 
C x ,C y ,C z values according to the scheme recommended in Section 6,7 and 8. 
This, for a specified satellite position, required computation of C x ,C y , C z as 
16 point mean for 2 blocks, 9 point mean for 3 blocks, 4 point mean for 11 
blocks and center point value for the remaining 168 blocks, and the time was 
0.055 seconds. 

It was generally felt, before this investigation was taken up, that the 
computation of C x , C y , C z would not be of adequate accuracy if done only at the 
center point of a 15° x 15° mean gravity anomaly block. The usual safeguard 
would than be to take a 4 point mean for all the 184 blocks. As we realize now, 
this would still give inaccurate results for the nearest few blocks, and would 
take 0.155 - 0.055 = 0. 1 seconds additional time for each satellite position. If 
we needed the computation at intervals of one minute for some orbital analysis 
studies, the extra CPU time for each 24 hour simulated orbit would be about 


28 



2. 5 minutes 


Let us consider the extreme case of the requirement as per Tables 
5 and 9, of computation of 16 point mean for 2 blocks, 9 point mean for 4 
blocks, 4 point mean for 21 blocks and center point value for 157 blocks. 
The saving in CPU execution time as compared to 4 point mean for all 
blocks would then still be 0. 09 seconds for each satellite position. At the 
same time, it will ensure that no errors larger than about 2. 5% are being 
caused for any block, which is not the case if 4 point mean is being taken 
for all blocks. 

10. Summary and Conclusions 


The computation of partial derivatives of the disturbing force of the 
earth's gravity field with respect to individual gravity anomalies is required 
to be computed for several orbital and trajectory analysis studies of artifical 
earth satellites used for geodetic purposes. We have here considered the 
case of solution of equation of motion of a satellite only affected by the earth's 
gravitational force, which is described in terms of l^ 3 x 15° equal area mean 
gravity anomaly blocks, as referred to a defined reference surface. For 
convenience, we define the partial derivatives of the components X , in an inertial 
coordinate system A., of the disturbing force C* , C y , C t given by: 



where T is the disturbing potential of the earth's gravity field and Agt is an 
individual mean gravity anomaly over one of the 184 15° x 15° equal area blocks. 
(For details see Section 3). 

We have examined in this report the numerical evaluation of C*,C y , C z . 

The principle criterion has been the spherical distance i|r * of the center of the 
block from the satellite subpoint. It has been found that for the blocks nearest 
to the satellite subpoint, the value of C x ,C y ,C z are in considerable error, even 
more than 50% for i|f* <10° (Tables 2 to 4) X, if computed only for the center of the 
block. This error reduces to about 2. 5% for i|r* exceeding 30° - 40° , depending 
on the height of the satellite (Tables 6 to 8), as compared to 16 point mean, 
described below. 
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For these blocks near the satellite subpoint, the value of C x ,C y ,C z can 
be computed within a tolerance of 1 to 3%, if the value is computed at several 
points in the block and then meaned. It has been found that these subdivisions 
of blocks should by symmetric with respect to latitude and longitude, thus 
giving 4 , 9 or 16 subdivisions and thus 4 point, 9 point and 16 point means 
for consideration . For I* <5?, 9 point mean is in error by 1 to 4% as 
compared to the 16 point mean. It therefore appears advisable to take 16 
point mean for a few of the nearest blocks. 

Four point means compare with 16 point means within about 2% for 
i|f* >15° to 20°, and for blocks nearer to the satellite subpoint, 9 point means 
require to be taken. The actual limiting value of \|i* for taking 16/9 point 
means may be chosen, in View of above figures, after considering the number 
of blocks which are likely to occur for i|»* from 5° to 20° at various latitudes. 
From Figure 1 and Table 5, the limiting values are found to be 10° for 16 point 
mean and 20° for 9 point mean. The maximum number of blocks for \|r*< 10° 
is two, which would have increased to four, if this limit was exceeded. Similary, 
the maximum number of blocks for t|r* <20° is six (inclusive of blocks for 16 
point mean) which would have increased to nine if this limiting value of f* was 
to exceed 20°. 

The computation of 16 point mean for '|/*<10° and 9 point mean for 
10° ^ i|f* •s 20° appears to hold for satellite heights from 400 km to 1600 km above 
the earth. For higher satellites, 16 point means may be dispensed with and 
9 point means may be taken for i|f* <20°. 

The limiting values of ijr* for taking 4 point means shows greater variation 
with the height of satellite, as discussed in Section 7. This has been found to 
be 20° '< ^ 30° for satellite height exceeding 1600 km, 20° < f ^ 35° for height 

from 800 - 1600 km, and 20° < tfr* ^ 45° for lower satellites. The total number 
of blocks corresponding to these limits would be 12 to 13, 12 to 18, 24 to 27 
respectively. These numbers include 4 to 6 blocks for which 9/16 point means 
would be taken. 

The values of C x , C y , C z do not decrease uniformly with increase in i|r*, 
and, in particular, it is not possible to dispense with the computation of values 
for those mean gravity anomaly blocks, which are near the antipode, the point 
at a spherical distance of \|t* = 180° from the satellite subpoint. (For details 
see Section 9). 

Finally, the actual CPU execution time was checked on IBM System/370 
Model 165 Computer for the 16/9/4 point mean and center point computation of 
C x ,C y ,Cz» as per limiting values of i|r* discussed above. It was found to have 
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taken 0. 1 seconds less for a single satellite position, as compared to the 4 
point mean for all the blocks.. iThe saving of time, for example, in the 
generation of satellite orbits over extended periods by numerical integration 
approach, would be noticeable. Further, the 4 point mean would have caused 
errors for the nearest blocks of 2 to 20%, while the errors in the scheme 
suggested now would be less than about 2. 5% for all blocks. 
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